Dynamics around supermassive black holes 



By ALESSIA GUALANDRIS and DAVID MERRITT 

Department of Physics and Center for Computational Relativity and Cravitation, Rochester 
Institute of Technology, Rochester, USA 

The dynamics of galactic nuclei reflects the presence of supermassive black holes (SBHs) in many 
ways. Single SBHs act as sinks, destroying a mass in stars equal to their own mass in roughly 
one relaxation time and forcing nuclei to expand. Formation of binary SBHs displaces a mass in 
stars roughly equal to the binary mass, creating low-density cores and ejecting hyper-velocity 
stars. Gravitational radiation recoil can eject coalescing binary SBHs from nuclei, resulting in 
offset SBHs and lopsided cores. We review recent work on these mechanisms and discuss the 
observable consequences. 



1. Characteristic scales 

Supermassive black holes (SBHs) are ubiquitous components of bright galaxies and 
many have been present with roughly their current masses (~ IO^AIq) since very early 
times, as soon as '--^ 1 Gyr after the Big Bang (Fan et al. 2003; Marconi et al. 2004). 
A SBH strongly influences the motion of stars within a distance r/j, the gravitational 
influence radius, where 

GM, 

rh^ — 1-1 

M, is the SBH mass and Uc is the stellar (Id) velocity dispersion in the core. Using the 
tight empirical correlation between M, and Uc'. 

( ^ ^ (1.66 ±0.24) ( , a = 4.86 ±0.4 (1.2) 

ViOSMq/ ^ ^V200 kms-V 

(Ferrarese & Ford 2005), this can be written 

2.86 . xO.59 



Th « 18 pc ^ « 13 pc • (1-3) 

VSOOkms^V KIO'^MqJ ^ ' 

While the velocities of stars must increase - by definition - inside r^, this radius is 
not necessarily associated with any other observational marker. Such is the case at the 
Galactic center, for instance, where the stellar density exhibits no obvious feature at 
Th ~ i pc. However the most luminous elliptical galaxies always have cores, regions 
near the center where the stellar density is relatively low. Core radii are of order rh in 
these galaxies, and the stellar mass that was (apparently) removed in creating the core 
is of order M,. These facts suggest a connection between the cores and the SBHs, and 
this idea has motivated much recent work, reviewed here, on binary SBHs and on the 
consequences of displacing SBHs temporarily or permanently from their central locations 
in galaxies. 

An important time scale associated with galactic nuclei (not just those containing 
SBHs) is the relaxation time, defined as the time for gravitational encounters between 
stars to establish a locally Maxwellian velocity distribution. The nuclear relaxation time 
is (Spitzer 1987) 
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Figure 1. Relaxation time, measured at the SBH influence radius, in a sample of early-type 
galaxies (Cote et al. 2004), vs. the central stellar velocity dispersion. Filled symbols are galaxies 
in which the SBH's influence radius is resolved; the star is the Milky Way bulge. (From Merritt, 
Mikkola & Szell 2007). 



1.2 X 10^° yr 



In A 



(1.46) 



lOOkms^V VlO^^/opc-^ 

with pc the nuclear density and hi A the Coulomb logarithm. Figure [1] shows estimates 
of Tr measured at Vh, in a sample of early- type galaxies, assuming = IM©. A least- 
squares fit to the points (shown as the dashed line in the figure) gives 

7.47 / 71 -r \ 1-54 



Tn(Tn) 



2.5 X 10^^ yr 



9.6 X 10^^ yr ' 



(1.5) 



.200 kms"V " ViO^Mq, 

"CoUisional" nuclei can be defined as those with T]i{rh) < 10 Gyr; figure [T] shows that 
such nuclei are uniquely associated with galaxies that are relatively faint, as faint as 
or fainter than the Milky Way bulge, which has Tjf{rf^) « 4 x 10^" yr. Furthermore, 
relaxation-driven changes in the stellar distribution around a SBH are generally confined 
to radii < 10~^rh, making them all but unobservable in galaxies beyond the Local Group 
(T. Alexander, these proceedings). But the relaxation time also fixes the rate of gravita- 
tional scattering of stars into the central "sink" - either a single or a binary SBH - and 
this fact has important consequences for nuclear evolution in low-luminosity galaxies, as 
discussed below. 



2. Core structure 

The 'core' of a galaxy can loosely be defined as the region near the center where the 
density of starlight drops significantly below what is expected based on an inward ex- 
trapolation of the overall luminosity profile. At large radii, the surface brightness profiles 
of early-type galaxies are well fit by the Sersic (1968) model, 

\nI{R) = In/e - b{n) [(i?/i?e)^^" ~ ll . (2.1) 
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The quantity b is normally chosen such that Re is the projected radius containing one- 
half of the total light. The shape of the profile is then determined by n; n = 4 is 



the de Vaucouleurs (1948) model, which is a good representation of bright elliptical (E) 
galaxies (Kormendy & Djorgovski 1989), while n = 1 is the exponential model, which 
approximates the luminosity profiles of dwarf elliptical (dE) galaxies (Binggeli, Sandage 
& Tarenghi 1984). An alternative way to write (|2.ip is 

din/ b f rV"" , , 

(2.2) 



d In R n \R, 

i.e. the logarithmic slope varies as a power of the projected radius. While there is no 
consensus on why the Sersic model is such a good representation of stellar spheroids, a 
possible hint comes from the dark-matter halos produced in A^-body simulations of hier- 
archical structure formation, which are also well described by (|2.2[) (Navarro et al. 2004), 
suggesting that Sersic's model applies generally to systems that form via dissipationless 
clustering (Merritt et al. 2005). 

Sersic's model is known to accurately reproduce the luminosity profiles in some galaxies 
over at least three decades in radius (e.g. Graham et al. 2003), but deviations often 
appear near the center. Galaxies fainter than absolute magnitude Mb « — 19 tend to 
have higher central surface brightness than predicted by Sersic's model; the structure 
of the central excess is typically unresolved but its properties are consistent with those 
of a compact, intermediate-age star cluster (Carollo, Stiavelli & Mack 1998; Cote et al. 
2006; Balcells et al. 2007). Galaxies brighter than Mb ~ —20 have long been known 
to exhibit central deficits (e.g. Kormendy 1985a); these have traditionally been called 
simply "cores," perhaps because a flat central density profile was considered a priori 
most natural (Tremaine 1997). 

For about two decades, it was widely believed that dE galaxies were distinct objects 
from the more luminous E galaxies. The dividing line between the two classes was put at 
absolute magnitude AIb ~ —18, based partly on the presence of cores in bright galaxies, 
and also on the relation between total luminosity and mean surface brightness (Kormendy 



1985b). This view was challenged by Jerjen & Binggeli (1997) and in a compelling series 



of papers, A. Graham and collaborators showed that - aside from the cores ~ early-type 
galaxies display a remarkable continuity of structural properties, from Mb ~ —13 to 
Mb « -22 (Graham & Guzman 2003; Graham et al. 2003; Trujillo et al. 2004). 

The connection between nuclear star clusters and SBHs, if any, is unclear; in fact it 
has been suggested that the two are mutually exclusive (Ferrarese et al. 2006; Wehner 
& Harris 2006), although counter-examples to this rule probably exist, e.g. NGC 3384 
which contains a nuclear cluster (Ravindranath et al. 2001) and may contain a SBH 
(Gebhardt et al. 2003). 

Here we focus on the cores. The cores extend outward to a break radius rf, that is 
roughly a few times r/j, or from 0.01 to ~ 0.05 times R^. A more robust way of 
quantifying the cores is in terms of their mass (i.e. light): the "mass deficit" (Milosavljevic 
et al. 2002) is defined as the difference in integrated mass between the observed density 
profile p(r) and an inward extrapolation of the outer profile, Poutif), typically modelled 
as a Sersic profile (figure [2]): 

Mdef = 4^ / [pout{r) - p{r)] r^dr. (2.3) 

JQ 

Figure [2] shows mass deficits for a sample of "core" galaxies, expressed in units of the 
SBH mass. There is a clear peak at Mdcf ~ IM,, although some galaxies have much 
larger cores. 
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Figure 2. Left: Surface brightness profile in the R band of NGC 3348, a "core" galaxy. The 
dashed line is the best-fitting Sersic model; the observed profile (points, and solid line) falls below 
this inside of a break radius rt « 0".35. (From Graham 2004.) Right: Histogram of observed mass 
deficits for the sample of core galaxies in Graham (2004) and Ferrarese et at. (2006). (Adapted 
from Merritt 2006a.) 



The fact that core and SBH masses are often so similar suggests a connection between 
the two. Ejection of stars by binary SBHs during galaxy mergers is a natural model 
(Begelman, Blandford & Rees 1980); the non-existence of cores in fainter galaxies could 
then be due to regeneration of a steeper density profile by star formation (e.g. McLaughlin 
et. at 2006) or by dynamical evolution associated with the (relatively) short relaxation 
times in faint galaxies (e.g. Merritt & Szell 2006). However the largest cores are difficult 
to explain via the binary model (Milosavljevic & Merritt 2001). 



3. Massive binaries 

A typical mass ratio for galaxy mergers in the local Universe is ^ 10 : 1 (e.g. Sesana 
et al. 2004) . To a good approximation, the initial approach of the two SBHs can therefore 
be modelled by assuming that the galaxy hosting the smaller BH spirals inward under the 
influence of dynamical friction from the fixed distribution of stars in the larger galaxy. 
Modelling both galaxies as singular isothermal spheres [p ~ r~^) and assuming that the 
smaller galaxy spirals in on a circular orbit, its tidally-truncated mass is m2 ~ a\r jlGa 
where ai and a are the velocity dispersion of the small and large galaxy respectively 
(Merritt 1984). Chandrasekhar's (1943) formula then gives for the orbital decay rate and 
infall time 

= -0.30 In A, t^nfaii « 3.3^^ 3.1 

dt ar (Tj 

where In A has been set to 2. Using (|1.2p to relate ct and (T2 to the respective SBH masses 
Ml and M2, this becomes 

Im/aH W 3.3-^ I — I , (3.2) 

i.e. tin fall exceeds the crossing time of the larger galaxy by a factor ^ 9^"'®, q = M2/M1 ^ 
1. Thus for mass ratios q > 10~^, infall requires less than ~ lO^Tcr ~ 10^" yr. This mass 
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= ~ -.jr—^Th (3.3a) 



ratio is roughly the ratio between the masses of the largest (~ lO^'^M©) and smallest 
(~ IO^-^Mq) known SBHs and so it is reasonable to assume that galaxy mergers will 
almost always lead to formation of a binary SBH in a time less than 10 Gyr. This 
conclusion is strengthened if the effects of gas are taken into account (e.g. Mayer et al. 
2007). 

Equation (|3.ip begins to break down when the two SBHs approach more closely than 
~ r/j, the influence radius of the larger hole, since the orbital energy of M2 is absorbed 
by the stars, lowering their density and reducing the frictional force. In spite of this 
slowdown, A^-body integrations (Merritt & Cruz 2001; Merritt & Milosavljevic 2001; 
Makino & Funato 2004; Berczik et al. 2005) show that the separation between the two 
SBHs continues to drop rapidly until the binary semi-major axis \s a k, Uh, where 

4^ ^ 4(l + g)2 

and n = M1M2/ {Ml + M2) is the binary reduced mass. At this separation - the "hard 
binary" separation - the binary's binding energy per unit mass is ~ and it ejects stars 
that pass within a distance ^ a with velocities large enough to remove them from the 
nucleus (Mikkola & Valtonen 1992; Quinlan 1996). 

What happens next depends on the density and geometry of the nucleus. In a spherical 
or axisymmetric galaxy, the mass in stars on orbits that intersect the binary is small, 
< M1-I-M2, and the binary rapidly ejects these stars; no stars then remain to interact with 
the binary and its evolution stalls (figure[3]). In non-axisymmetric (e.g. triaxial) nuclei, the 
mass in stars on centrophilic orbits can be much larger, allowing the binary to continue 
shrinking past a^- And in coUisional nuclei of any geometry, gravitational scattering of 
stars can repopulate depleted orbits. These different possbilities are discussed in more 
detail below. 

If the binary does stall at a « a/j, it will have given up an energy 

^^^GM^^GM^ (3.4.) 
2rh 2ah 

«_iA/2(j2 + 2(Mi + M2)a2 (3.46) 

« 2(Mi + Af2)cr^ (3.4c) 

to the stars in the nucleus, i.e., the energy transferred from the binary to the stars is 
roughly proportional to the combined mass of the two SBHs. The reason for this counter- 
intuitive result is the Uh oc M2 dependence of the stalling radius p.3p : smaller infalling 
BHs form tighter binaries. Detailed A^-body simulations (Merritt 2006a) verify that the 
mass deficit generated by the binary in evolving from to at is a weak function of 
the mass ratio, 

q 



Mdof,h « 0.4(Mi + M2) [-^j (3.5) 

for 0.025 < (7 < 0.5. A mass deficit of ~ 0.5M, is still a factor ^ 2 too small to explain the 
observed peak in the Mdd/M, histogram (figure [2]). On the other hand, bright elliptical 
galaxies have probably undergone numerous mergers, and the proportionality between 
Mdef and Mi + M2 (rather than, say, M2) implies that the mass deficit following J\f 
mergers is ^ 0.5J\f times the accumulated BH mass. Mass deficits in the range 0.5 < 
Mdef/M, < 1.5 therefore imply 1 ^ A/" < 3 mergers, consistent with the number of 
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Figure 3. Early evolution of binary SBHs in A'^-body galaxies, for different values of the binary 
mass ratio, M2/M2 — 0.5, 0.25, 0.1, 0.05, 0.025 (left to right). The upper horizontal line indicates 
rh, the influence radius of the more massive hole. Lower horizontal lines show (|3.3p . the 
"hard- binary" separation. The evolution of the binary slows drastically when a ~ a^; in a real 
(spherical) galaxy with much larger A'^, evolution would stall at this separation. The smaller the 
infalling BH[, the farther it spirals in before stalling. (Adapted from Merritt 2006a.) 



major mergers expected for bright galaxies since the epooch at which most of the gas 
was depleted (e.g. Haehnelt & KaufFmann 2002). Hierarchical growth of cores tends to 
saturate after a few mergers however making it difficult to explain mass deficits greater 
than ^ 2M, in this way. An effective way to enlarge cores still more is to kick the SBH 
out, at least temporarily, as discussed in §4. 

The first convincing evidence for a true, binary SBH was recently presented by Rodriguez et al. (2006)1 
who discovered two compact, flat-spectrum AGN at the center of a single elliptical galaxy, 
with a projected separation of 7 pc. This is consistent with the radius at formation of 
a ~ 1O''-^M0 binary (|1.3p . or the stalling radius for a binary of at least ~ lO^ '^M© p.3p . 
Rodriguez et al. estimate a binary mass of ~ IO^Mq but with considerable uncertainty. 
All other examples of "binary" SBHs in single galaxies have separations ^ (Komossa 
2006). 

Even in a spherical galaxy, the stalling that occurs at a w a/i can be avoided if stars 
continue to be scattered onto orbits that intersect the binary (Valtonen 1996; Yu 2002; 
Milosavljevic & Merritt 2003). Such "coUisional loss-cone repopulation" requires that 
the two-body (star-star) relaxation time at r rh be less then ~ 10^" yr; according to 
(|1.5|) . this is the case in galaxies with Af, < 10^ Mq, i.e. at the extreme low end of the 
SBH mass distribution. CoUisional loss cone repopulation is therefore irrelevant to the 
luminous galaxies that are observed to have cores but may be important in the mass 
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t t 
Figure 4. Short-term (a) and long-term (b) evolution of a massive binary in a series of A'^-body 
integrations. Vertical axis is the inverse semi-major axis (i.e. energy) of the binary, computed 
by averaging several independent A'^-body runs; different curves correspond to different values of 
A'', the number of "star" particles. The evolution of the binary is independent of until a Pi an 
(horizontal line) ; thereafter the evolution rate is limited by how quickly stars are scattered onto 
orbits that intersect the binary, and decreases with increasing N . (From Merritt, Mikkola & 
SzeU 2007.) 

range (M, < lO^M©) of most interest to space-based gravitational wave interferometers 
like LISA (Hughes 2006). 

TV-body simulation would seem well suited to this problem (e.g. Governato, Colpi & 
Maraschi 1994; Makino 1997; Milosavljevic & Merritt 2001; Makino & Funato 2004). 
The difficulty is noise - or more precisely, getting the level of noise just right. In a real 
galaxy there is a clear separation of time scales between an orbital period, and the time 
for stars to be scattered onto depleted orbits: the first is typically much shorter than 
the second which means that orbits intersecting the binary will remain empty for many 
periods before a new star is scattered in. This is called the "empty loss cone" regime 
and it implies that supply of stars to the binary will take place diffusively. In an TV- 
body simulation, however, N is much smaller than its value in real galaxies and orbits 
are repopulated too quickly. This is the reason that binary evolution rates in TV-body 
simulations typically scale as TV^",a < 1 rather than the ~ T^^ ~ TV~^ dependence 
expected if stars diffused gradually into an empty loss cone (Merritt & Milosavljevic 
2003). Figure m provides an illustration: early evolution of the binary, until a « a^, is 
TV- independent; formation of a hard binary then depletes the loss cone and continued 
hardening occurs at a rate that is a decreasing function of TV, though less steep than 
TV-i. 

An alternative approach is based on the Fokker-Planck equation. Both single and 
binary SBHs can be modelled as "sinks" located at the centers of galaxies (Yu 2002). 
The main differences are the larger physical extent of the binary (~ G{Mi + M2)/a'^ vs. 
GMt/c^) and the fact that the binary gives stars a finite kick rather than disrupting or 
consuming them completely. However the diffusion rate of stars into a central sink varies 
only logarithmically with the size of the sink (Lightman & Shapiro 1978), and a hard 
binary ejects most stars well out of the core with y 2> cr, so the analogy is fairly good. 




Figure 5. Evolution of a binary SBH in a collisional nucleus, based on a Fokker-Planck model 
that allows for evolution of the stellar distribution (Merritt, Mikkola & Szell 2007). Left: 
Probability of finding the binary in a unit interval of In a. From left to right, curves are for 
Ml +M2 = (0.1,1,10,100) X IO^Mq. Solid (dashed) curves are for M2/M1 = 1(0.1). Open 
circles indicate when the rate of energy loss to stars equals the loss rate to gravitational waves; 
filled circles correspond to an elapsed time since a = of 10^° yr. For the two smallest values 
of A/,, the latter time occurs off the graph to the right. Right: Total time for a binary to evolve 
from a = ah to gravitational wave coalescence as a function of binary mass. The thick (black) 
curve is for M2/M1 — 1 and the thin (blue) curve is for M2/M1 = 0.1. Dotted curves show the 
time spent in the gravitational radiation regime only. 



The Fokker-Planck equation describing nuclei with sinks is (Bahcall & Wolf 1977) 

Here N{E,t) = ATTp{E)f{E,t) is the distribution of stellar energies, f{E,t) is the phase 
space density and p{E) is a phase-space volume element. The first term on the RHS of 
(|3.6p describes the response of / to the flux Fe of stars in energy space due to encounters. 
The second term, —J^, is the flux of stars into the sink, which is dominated by scattering 
in angular momentum (Frank & Rees 1976). A proper treatment of the latter term 
requires a 2d (energy, angular momentum) analysis, but a good approximation to T can 
be derived by assuming that the distribution of stars has reached a quasi steady-state 
near the loss cone boundary in phase space (Cohn & Kulsrud 1979). If the sink is a 
binary SBH, a second equation is needed that relates the flux of stars into the loss cone 
to the rate of change of the binary's semi-major axis: 



d fl\ 2(C) 



dt \a J a(Mi -I- M2 



T{E,t)dE (3.7) 



with (C) ~ 1.25 a dimensionless mean energy change for stars that interact with the 
binary. 

Both terms on the RHS of (|3.6p imply changes in a time ~ Tr. The first term on its 
own implies evolution toward the Bahcall- Wolf (1976) "zero flux" solution, p ~ r^'^l'^. 
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The second term implies that a mass of order A/, will be scattered into the sink in 
a time of T[i{rh). When the sink is a binary SBH, the binary responds by ejecting the 
incoming stars and shrinking, according to (|3.7p . As a result, changes in the structure of 
the nucleus on a relaxation time scale (e.g. growth of a core) are directly connected to 
changes in the binary semi-major axis. 

Numerical solutions to p.6p , p.7p (including also the effects of a changing gravitational 



potential) have been presented by Merritt et al. (2007) The solutions are well fit by 



, /■ah\ -D z t , , 

where t is defined as the time since the binary first became hard (a = a^), and the 
coefhcients A « 0.016, B « 0.08 depend weakly on the binary mass ratio. Including the 
effect of energy lost to gravitational radiation: 



d /1\ d /1\ d (\ 



dt \ a I dt \ a / dt 



(3.9) 



GR 



allows one to compute the time to full coalescence, Tcoai. Figure [5] shows Tcoai (right 
panel), and the time spent by the binary in unit intervals of In a prior to coalescence (left 
panel), as functions of binary mass. The time to coalescence is well fit by 



^-log.l^^V (3.10c) 



with 



Y ^Ci+C2X + C3X', (3.10a) 

-^coal 
lO^Vy 

lO^Me 

M2/M1 = 1 : Ci = -0.372, C2 = 1.384, C3 = -0.025 (3.11a) 
M2/M1 = 0.1 : Ci = -0.478, C2 = 1.357, C3 = -0.041. (3.11&) 

Based on the figure, binary SBHs would be able to coalesce via interaction with stars 
alone in galaxies with M, < 2 x lO^^M©. For A/1 + M2 > lO^M©, evolution for 10 Gyr only 
brings the binary separation slightly below ah', in such galaxies the most likely separation 
to find a massive binary (in the absence of other sources of energy loss) would be near 
ah- 

The core continues to grow as the binary shrinks, but the mass deficit is not related in 
a simple way to the mass in stars "ejected" by the binary (e.g. Quinlan 1996). Rather it 
results from a competition between loss of stars to the binary, represented by ~J-{E,t), 
and the change in N{E, t) due to diffusion of stars in energy, represented by —dFE/dE. 
As the mass deficit increases, so do gradients in /, which increases the fiux of stars toward 
the center and counteracts the drop in density. In principle the two terms could balance, 
but at some distance from the center the relaxation time is so long that local Fe (E) must 
drop below the integrated loss term !F{E)dE ~ stars can not diffuse in fast enough to 
replace those being lost to the binary and the density drops. The Fokker-Planck solutions 
show that the mass deficit increases with binary binding energy as 

Afdef.c ~ 1.7 {Ml + M2) logio [ah/a) (3.12) 



again with a weak dependence on M2/M1. The mass deficit at the onset of the gravita- 
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tional radiation regime is found to be 

Mdcf,c « (4.5, 3.5, 2.6, 1.6)(Mi + Ms) {M2/M1 = 1) (3.13a) 
« (3.4,2.6,1.7,0.9)(Afi + M2) (M2/A/1 = 0.1) (3.13&) 

for M1 + M2 = (10^, 10*^, 10^, 10^)Mq. These values should be added to the mass deficits 
(|3.5p generated during formation of the binary when predicting core sizes in real galaxies. 

Are such mass deficits observed? Only a handful of galaxies in the relevant mass range 
(Mgai < 10^°Mq) are near enough that their cores could be resolved even if present; of 
these, neither the Milky Way nor M32 exhibit cores. Also, as noted above, many low- 
luminosity spheroids have compact central excesses rather than cores. These facts do 
not rule out the past existence of massive binaries in these galaxies however. (1) Binary 
evolution might have been driven more by gas dynamical torques than by ejection of 
stars; gas content during the most recent major merger is believed to be a steep inverse 
function of galaxy luminosity (Kauffmann & Haehnelt 2000). (2) Star formation can 
create a dense core after the two SBHs have coalesced (Mihos & Hernquist 1994). (3) A 
two-body relaxation time short enough to bring the two SBHs together would also allow 
a Bahcall-Wolf cusp to be regenerated in a comparable time after the two SBHs combine, 
tending to erase the core (Merritt & Szell 2006). 

From the point of view of physicists hoping to detect gravitational waves, it is disap- 
pointing that this model only guarantees coalescence at the extreme low end of the SBH 
mass distribution. (Astronomers hoping to detect binary SBHs may take the opposite 
point of view.) Fortunately, there is no dearth of ideas for overcoming the "final parsec 
problem" and allowing binary SBHs to merge efficiently, even in massive galaxies: 

Non-£Lxisymmetric geometries. Real galaxies are not spherical nor even axisym- 
metric; parsec-scale bars are relatively common and departures from axisymmetry are 
often invoked to enhance fueling of AGN (e.g. Shlosman et al. 1990). Orbits in a triaxial 
nucleus can be "centrophilic," passing arbitrarily close to the center after a sufficiently 
long time (Poon & Merritt 2001, 2004). This implies feeding rates for a central binary 
that can approach the "full loss cone" rate in spherical geometries, or 

if the fraction of centrophilic orbits is large (Merritt & Poon 2004). While is 
impossible to know in any particular galaxy, even small values imply much larger feeding 
rates than in a diffusively-repopulated loss cone. Figure [S] shows results from A'^-body 
simulations that support this idea. 

Secondary slingshot. Stars ejected by a massive binary can interact with it again 
if they return to the nucleus on nearly-radial orbits. The total energy extracted from 
the binary via this "secondary slingshot" will be the sum of the discret energy changes 



during the interactions. Milosavljevic & Merritt (2003) showed that a mass Mi, of stars 



initially in the binary's loss cone causes the binary to evolve as 

1 1 4 / t\ 2ucr2 , ^ , , 

-«- + -ln 1 + - , ^0= ,/i.r,. PiE) 3.15 
a ah rh \ hj Mi,[AE} 

in the absence of diffusive loss cone repopulation, where (A£') is the specific energy change 
after one interaction with the binary, E is the initial energy and P{E) is the orbital period. 



The secondary slingshot runs its course after a few orbital periods. Sesana et al. (2007) 
sharpened this analysis by carrying out detailed three-body scattering experiments and 
recording the precise changes in energy of stars as they underwent repeated interactions 
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Figure 6. Efficient merger of binary SBHs in barred galaxies. Plots are based on N-body 
simulations (no gas) of equal-mass binaries at the centers of galaxy models, with and without 
rotation, (a) Spherical models. The binary hardening rate declines with increasing A^, as in 
figure [4l implying that the evolution would stall in the large- A'^ limit, (b) Binary evolution in 
a flattened, rotating version of the same galaxy model. At t 10, the rotating model forms a 
triaxial bar. Binary hardening rates in this model are essentially independent of A'^, indicating 
that the supply of stars to the binary is not limited by coUisional loss-cone refilling as in the 
spherical models. This is currently the only simulation that follows two SBHs from kiloparsec 
to sub-parsec separations and that can be robustly scaled to real galaxies. (From Berczik et at. 
2006.) 



with the binary. They inferred modest (~ x2) changes in 1/a due to the secondary shng- 
shot, but their assumption of a p r^^ density profile around the binary was probably 
over-optimistic; such steep density profiles are never observed and even if present initially 
would be rapidly destroyed when the binary first formed. 

Bound subsystems. As noted above, recent observational studies have greatly in- 
creased the number of galaxies believed to harbor compact nuclear star clusters; inferred 
masses for the clusters are comparable with the mass that would normally be associated 
with a SBH. It is not yet clear whether these subsystems co-exist with SBHs, but if 
they do, they could provide an extra source of stars to interact with a massive binary. 
Zier (2006, 2007) explored this idea, assuming a steeply rising density profile around the 
binary, p oc r"^, at the time that its separation first reached ah- Zier concluded that a 
cluster having total mass ~ Mi -I-M2, distributed as a steep power law, 7 > 2.5, could ex- 
tract enough energy from the binary to allow gravitational wave coalescence in less than 



10 Gyr. iV-body tests of this hypothesis are sorely needed; as in the Sesana et al. (2007) 



study, Zier's approach did not allow him to self-consistently follow the effect of formation 
of the binary on the surrounding mass distribution. 

Masive perturbers. In a nucleus containing a spectrum of masses, the gravitational 
scattering rate is proportional to 

J n{m)rn^dm 



J n(m)mdm 



(3.16) 



(e.g. Merritt 2004). Perets et al. (2007) argued that "massive perturbers" near the center 
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of the Milky Way - massive stars, star clusters, giant molecular clouds - are sufficiently 
numerous to dominate m, implying potentially much higher rates of gravitational scatter- 



ing into a central sink than in the case of solar-mass perturbers. Perets & Alexander (2007) 



extended this argument to galaxies in general, emphasizing in particular the early stages 
following a galactic merger, and concluded that coUisional loss cone repopulation would 
be sufficient to guarantee coalescence of binary SBHs in less than 10 Gyr for all but 
the most massive binaries. As in the studies of Sesana et al. (2007)j and Zier (2007) 



Perets & Alexander (2007) optimistically assumed a steep (p oc r ) density profile 



around the binary, in spite of iV-body studies showing rapid destruction of the cusps. 
Their arguments for massive perturbers in giant E galaxies are also rather speculative. 

Multiple SBHs. An extreme case of a "massive perturber" is a third SBH, which 
might scatter stars into a central binary (Zhao et al. 2002), or perturb the binary di- 
rectly, driving the two SBHs into an eccentric orbit and shortening the time scale for 
gravitational wave losses (Valtonen et al. 1994; Makino & Ebisuzaki 1994; Blaes et al. 
2002; Volonteri et al. 2003; Iwazawa et al. 2006; Hoffman & Loeb 2007). The hkelihood 
of multiple-SBH systems forming is probably highest in the brightest E galaxies since 
massive binaries are most likely to stall (low stellar density, little gas) and since large 
galaxies experience the most frequent mergers. Here again, more iV-body simulations, 
including post-Newtonian terms, are needed; among other dynamical effects that could 
then be self-consistently included are changes in core structure, and BH-core oscillations 
like those described in the next section. 

Gas. The same galaxy mergers that create binary SBHs can also drive gas into the 
nucleus, and there is abundant observational evidence for cold (e.g. Jackson et al. 1993; 
Gallimore et al. 2001; Greenhill et al. 2003) and hot (e.g. Baganoff et al. 2003) gas near the 
centers of at least some galaxies. Dense concentrations of gas can substantially accelerate 
the evolution of a massive binary by increasing the drag on the individual BHs (Escala 
et al. 2004, 2005; Dotti et al. 2007). The plausibility of such dense accumulations of gas, 
with mass comparable to the mass of the SBHs, is unclear however (e.g. Sakamoto et al. 
1999; Christopher et al. 2005). Large-scale galaxy merger simulations (Kazantzidis et al. 
2005; Mayer et al. 2007) show that the presence of gas leads to more rapid formation 
of the massive binary, but these simulations still lack the resolution to follow the binary 
past a « r/j and so have nothing relevant to say (yet) concerning the final parsec problem. 



As this summary indicates, many possible solutions to the "final parsec problem" exist, 
but none is guaranteed to be effective in all or even most galaxies. The safest bet is that 
both coalesced and uncoalesced binary SBHs exist, but with what relative frequency is 
still anyone's guess. 



4. SBH/IBH binaries 

Secure dynamical evidence exists for SBHs in the mass range 10^'^ < M,/Mq < 10^'^ 
(Ferrarese & Ford 2005) and compelling arguments have been made for BHs with masses 
10^ < M,/Mq < 10 '' in active nuclei (Greene & Ho 2004). Binary mass ratios as extreme 
as 1000 : 1, and possibly greater, are therefore to be expected. This possibility has 
received most attention in the context of intermediate-mass black holes (IBHs) in the 
Milky Way, where they could form in dense star clusters like the Arches or Quintuplet 
before spiralling into the center and forming a tight binary with the ^ 3.5 x lif'MQ SBH 
(Portegies Zwart & McMillan 2002; Hansen & Milosavljevic 2003). 
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Figure 7. Left: A^-body simulations of the inspiral of IBHs into the center of the Milky Way. 
Solid lines show the separation between the two BHs and dashed lines are theoretical predictions 
that ignore loss-cone depletion or changes in the structure of the core. Smaller IBHs spiral 
in farther before "stalling"; the Mibh = 10^ Mq simulation ends before the stalling radius 
is reached. (From Baumgardt, Gualandris & Portegies Zwart 2006.) Right: Evolution beyond 
a = ah, based on the Fokker-Planck model of Merritt, Mikkola & Szell (2007). Dashed lines 
indicate when the evolution time due to gravitational radiation losses is less than 10 Gyr. 



The predicted hard-binary separation for a SBH/IBH pair is 

M. Mibh 



ah 



O.Smpc ^ 



10-^ 



3.5 X 10<^Mp 



q ■ 



M, 



(4.1) 



This separation - ^ 10^ AU - is comparable to the periastron distances of the famous 
"S" stars (Eckart et al. 2002; Ghez et al. 2005). Dynamical constraints on the existence of 
an IBH at this distance from the SBH are currently weak (Yu & Tremaine 2003; Hansen 
& Milosavljevic 2003; Reid & Brunthaler 2004). Figure [7^ shows A^-body simulations 
designed to mimic inspiral of IBHs into the Galactic center. The figure confirms the 
expected slowdown in the inspiral rate at a separation ~ a^. Figure [7)3 plots evolutionary 
tracks for the same three IBH masses as in the left panel, based on the Fokker-Planck 
model of Merritt et al. (2007) For Mibh ^ IQ^Mq, evolution of the binary is dominated 
by gravitational wave losses already at a = ah- 

The same inward flux of stars that allows the binary to shrink also implies an outward 
flux of stars ejected by the binary. The latter are a potential source of "hyper-velocity 
stars" (HVSs), stars moving in the halo with greater than Galactic escape velocity (Hills 
1988). The relation between the stellar ejection rate and the binary hardening rate, when 
a ^ tth, is given by (|3.7p after rewriting it as 



d (\ 



ard 



a—r 

dt 



2.5 

m7 



X flux; 



(4.2) 



here "flux" is the total mass in stars per unit time, from all energies, that are scattered 
into (and ejected by) the binary. Combining (|4.2[) with (|3.8p . the flux is 



5.0 



M. 

TR{rh) 



1 



5t 

TR{rh)_ 



-1-1/2 



(4.3a) 
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Figure 8. Left: Creation of a core by inspiral of IBHs into the Galactic center. The initial 
density profile is shown by the solid line. Green (dotted) line shows the core 10 Myr after 
the W^Mq IBH has merged with the SBH; almost no change occurs during this time. (From 
Baumgardt, Gualandris & Portegies Zwart 2006.) Right: Fokker-Planck model showing how the 
cusp regenerates due to two-body scattering, on Gyr timescales. (Adapted from Merritt, Mikkola 
& SzeU 2007.) 



where the second line uses values appropriate to the Galactic center. Relating the to- 
tal ejected flux to the number of HVSs that would be observed is not straighforward; 
for instance, only a fraction (< 10%) would be ejected with high enough velocity to 
still be moving faster than ~ 500 km s~^ after climbing through the Galactic potential 
(Gualandris et al. 2005; Baumgardt et al. 2006), and targeted searches for HVSs only 
detect certain stellar types so that knowledge of the stellar mass function is also required 
(Brown et al. 2006). 

Inspiral of the IBH creates a core of radius ~ 0.05 pc « 1" (figure [5^). Such a core 
might barely be detectable at the center of the Milky Way from star counts. There is 
no clear indication of a core (Schoedel et. al 2007), but if the inspiral occurred more 
than a few Gyr ago, star-star gravitational scattering would have gone some way toward 
"refilling" the region depleted by the binary (Merritt & Wang 2005; Merritt & SzeU 2006; 
figure [TJi). In this case, however, the ejected stars would almost all have moved beyond 
the range of HVS surveys by now. 

The angular distribution of the ejected stars has been proposed as a test for their origin; 
unlike other possible sources of HVSs, a SBH/IBH binary tends to eject stars parallel to 
the orbital plane or, if the orbit is eccentric, in a particular direction (Levin 2006; Sesana 
et al. 2006). In two A^-body simulations of IBH inspiral however (Baumgardt, Gualandris 
& Portegies Zwart 2006; Matsubayashi et al. 2007), the orientation of the binary began 
to change appreciably, in the manner of a random walk, after it became hard. This was 
due to "rotational Brownian motion" (Merritt 2002): torques from passing stars - the 
same stars that extract energy and angular momentum from the binary - also change 
the direction of the binary's orbital angular momentum vector. In one hardening time 



< 350MQyr 



-1 



(4.36) 
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\a/a\ of the binary, its orientation changes by 

■ VlO-3/ yiO^mJ 5 

(The eccentricity dependence in (j4.4&P is approximate; the numerical coefficient in this 
equation has only been confirmed by detailed scattering experiments for e = 0.) In both 
of the cited iV-body studies, the binary eccentricity evolved appreciably away from zero 
before the orientation changes became signficant. Rotational Brownian motion might 
not act quickly enough to randomize the orienation of a SBH/IBH binary in a time of 
~ 10® yr, the flight time from the Galactic center to the halo, unless perturbers more 
massive than Solar-mass stars are present near the binary however (Merritt 2002; Perets 
& Alexander 2007). 



(4.4a) 
(4.46) 



5. Kicks and cores 

After seeming to languish for several decades, the field of numerical relativity has re- 
cently experienced exciting progress. Following the breakthrough papers of Pretorius (2005)] 



Campanelli et al. (2006) and Baker et al. (2006a)| several groups have now successfully 



simulated the evolution of binary BHs all the way to coalescence. The final inspiral is 
driven by emission of gravitational waves, and in typical (asymmetric) inspirals, a net 
impulse - a "kick" - is imparted to the system due to anisotropic emission of the waves 
(Bekenstein 1973; Fitchett 1984; Favata et al. 2004). Early arguments that the magni- 
tude of the recoil velocity would be modest for non-spinning BHs were confirmed by the 
simulations, which found 14ick ^ 200 km s~^ in the absence of spins (Baker et al. 2006; 
Gonzalez et al. 2007a; Herrmann et al. 2007). The situation changed dramatically fol- 
lowing the first (Campanelli et al. 2007a) simulations of "generic" binaries, in which the 
individual BHs were spinning and tilted with respect to the orbital angular momentum 
vector. Kicks as large as ^ 2000 km s^^ have now been confirmed (Campanelli et al. 
2007b; Gonzalez et al. 2007b; Tichy & Marronetti 2007), and scaling arguments based 
on the post-Newtonian approximation suggest that the maximum kick velocity would 
probably increase to ~ 4000 km s^^ in the case of maximally-spinning holes (Campanelli 
et al. 2007b) . The most propitious configuration for the kicks appears to be an equal-mass 
binary in which the individual spin vectors are oppositely aligned and oriented parallel 
to the orbital plane. For unequal-mass binaries, the maximum kick is 

Knax ~ 6 X lO^km (5.1) 

(l + qP 

where q = M2/M1 ^ 1 is the binary mass ratio and maximal spins have been assumed 
(Campanelli et al. 2007c). Orienting the BHs with their spins perpendicular to the orbital 
angular momentum may seem odd (Bogdanovic, Reynolds & Miller 2007), but there is 
considerable evidence that SBH spins bear no relation to the orientations of the gas disks 
that surround them (e.g. Kinney et al. 2000; Gallimore et. al. 2006) and this is presumably 
even more true with respect to the directions of infalling BHs. Galaxy escape velocities 
are < 3000 km s^^ (Merritt et al. 2004), so gravitational wave recoil can in principle 
eject coalescing SBHs completely from galaxies. 

Detailed A^-body simulations show that the motion of a SBH that has been kicked 
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Figure 9. Left: Core oscillations in an A*'-body simulation of ejection of a SBH from the center 
of a galaxy; the kick velocity was 60% of the escape velocity. Contour plots show the stellar 
density at equally spaced times, spanning ~ 1/2 of the SBH's orbital period. Filled circles are 
the SBH and crosses indicate the location of the (projected) density maxima. (From Gualandris 
& Merritt 2007.) Right: Surface brightness contours of three "core" galaxies with double or offset 
nuclei, from Lauer et al. (2005). Top: NGC 4382; middle: NGC 507; bottom: NGC 1374. 

with enough velocity to eject it out of the core, but not fast enough to escape the galaxy 
entirely, exhibits three distinct phases (Gualandris & Merritt 2007): 

• Phase I: The SBH oscillates with decreasing amplitude, losing energy via dynamical 
friction each time it passes through the core. Chandrasekhar's theory accurately repro- 
duces the motion of the SBH in this regime for values 2 < InA < 3 of the Coulomb 
logarithm, if the gradually-decreasing core density is taken into account. 

• Phase II: When the amplitude of the motion has decayed to roughly the core radius, 
the SBH and core begin to exhibit oscillations about their common center of mass (fig- 
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Figure 10. Evolution of the SBH kinetic energy following a kick of 60% the central escape veloc- 
ity, in two TV-body simulations of a galaxy represented by A'^ stars, with A'^ = (2.5 x 10^, 2 x 10^). 
The mass of the SBH, and the total mass of the galaxy, are the same in the two simulations; 
all that varies is the mass of the "star" particles. The right-hand panel shows binned values of 
V^. Most of the elapsed time is spent in SBH/core oscillations like those illustrated in figure [O] 
Eventually, the SBH's kinetic energy decays to the Brownian value, shown as the horizontal 

1/2 

dashed lines in the right panel. The Brownian velocity scales as and so is smaller for larger 
A''. Scaled to a 3 x W^^Mq galaxy, the time to reach the Brownian regime would be ~ 10* yr. 
(Adapted from Gualandris & Merritt 2007.) 



ure [9]) . These oscillations decay exponentially (figure [T0| , but with a time constant that 
is 10 — 20 times longer than would be predicted by a naive application of the dynamical 
friction formula. 

• Phase III: Eventually the SBH's kinetic energy drops to an average value 

Im.K' « im.vj (5.2) 

i.e. to the kinetic energy of a single star. This is the regime of gravitational Brownian 
motion (Bahcall & Wolf 1976; Young 1977; Merritt et al. 2007). 

A natural definition of the "return time" of a kicked SBH is the time to reach the 
Brownian regime. Unless the kick is very close to the escape velocity, the return time is 
dominated by the time spent in "Phase II;" during this time, the SBH's energy decays 
roughly as 

£: « $(0) + $(re)e-(*-*'^)/^ (5.3) 

(Gualandris & Merritt 2007); tc is the time when the SBH re-enters the core whose radius 
is Tc. The damping time in the A^-body simulations, t, is 

« 1.2 X lOV ( ^ ' { ^^^1 ( , (5.4&) 

^ V 250km s-V VlO'^^oPC^V VIO"^/©/ 

with pc and CTc the core density and velocity dispersion respectively. The number of decay 
times required for the SBH's energy to reach the Brownian level is ~ ln(M,/TO*) « 20, 
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Figure 11. Left: Points show projected density profiles computed from A''-body models af- 
ter a kicked SBH has returned to the center, for three different values of the kick velocity 
(Vkick = (0.2,0.4,0.8) X Ksc)- Each set of points is compared with the best-fitting core-Sersic 
model (lines). The insert shows a zoom into the central region. Right: Mass deficits generated 
by kicked SBHs. The differents sets of points correspond to different galaxy models. (Adapted 
from Gualandris & Merritt 2007.) 



implying that a kicked SBH will remain significantly off-center for a long time, as long 
as 1 Gyr in a bright galaxy with a low-density core. 

In fact, asymmetric cores are rather common. These include off-center nuclei (Bin- 
gelli et al. 2000; Lauer et al. 2005); double nuclei (Lauer et al. 1996); and cores with 
a central minimum in the surface brightness (Lauer et al. 2002). Three examples, from 
Lauer et al. (2005)] are reproduced here on the right side of figure [9l all are luminous 
"core" galaxies, and each strikingly resembles at least one frame from the TV-body mon- 
tage on the left. The longevity of the "Phase 11" oscillations makes the kicks a plausible 
model for the observed asymmetries. This explanation is probably not appropriate for 
the famous double nucleus of M31, since M31 is not a "core" galaxy, and since one of 
the brightness peaks in M31 (the one associated with the SBH) lies essentially at the 
galaxy photocenter; Figure [9] suggests that an oscillating SBH would typically (though 
not always) be found on the opposite side of the galaxy from the point of peak brightness. 
The M31 double nucleus has been successfully modelled as a clump of stars on eccentric 
orbits which maintain their lopsidedness by virtue of moving deep within the Keplerian 
potential of the SBH (Tremaine 1995). 

The kicks are quite effective at inflating cores (Merritt et al. 2004; Boylan-Kinchin 
et al. 2004). FigurefTTl from' Gualandris fc Merritt (2007)] illustrates this: the mass deficit 
generated by the kick is approximately 

Mdof,k « 5Af. (Fkick/Ksc)'-^' . (5.5) 

Even modest kicks can generate cores substantially larger than those produced during 
the formation of a massive binary (~ IM,; eq. 13. 5p . Furthermore, this mechanism is 
potentially effective in even the most luminous galaxies, unlike the relaxation-driven 
model for core growth discussed above (eq. I3.12p which only applies to nuclei with short 
relaxation times. Gravitational radiation recoil is therefore a tenable explanation for 
the subset of luminous E galaxies with large mass deficits (figure [2]). An alternative 
explanation for the over-sized cores (Lauer et al. 2007) postulates that the SBHs in these 
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galaxies are "hypermassive," M, > lO^^M© and that the cores are a consequence of 
shngshot ejection by a massive binary. 

The kicks have a number of other potentially observable consequences, including spa- 
tially and/or kinematically offset AGN (Madau & Quataert 2004; Haehnelt et al. 2006; 
Merritt et al. 2006; Bonning et al. 2007) and distorted or wiggling radio jets (Gualandris 
& Merritt 2007). Many of these manifestations were first discussed by R. Kapoor in a 
remarkably prescient series of papers (Kapoor 1976, 1983a, b, 1985). 



6. Black-hole-driven expansion 

The growth of a core around a shrinking, binary SBH was discussed above: beyond 
a certain radius, the relaxation time becomes so long that the encounter-driven flux of 
stars toward the center cannot compensate for losses to the binary, forcing the density 
to drop. A similar process takes place around a single SBH (Shapiro 1977; Dokuchaev 
1989): stars coming too near are consumed, or disrupted, and the density drops. This 
effect is absent from the classical equilibrium models for stars around a BH (e.g. Bahcall 
& Wolf 1976; Cohn & Kulsrud 1978) since these solutions fix the phase-space density far 
from the BH, enforcing an inward flux of stars precisely large enough to replace the stars 
being consumed by the sink. In reality, the BH acts as a heat source, in much the same 
way that hard binary stars inject energy into a post-core-collapse globular cluster and 
cause it to re-expand. 

A simple model that produces self-similar expansion of a nucleus containing a SBH can 



be constructed by simply changing the outer boundary condition in the^Bahcall fc Wolf (1976) 



problem from /(O) = /o to /(G) = 0. One finds that the evolution after ~one relax- 
ation time can be described as p{r,t) = pc{t)p*{r), with p*{r) slightly steeper than the 
p ~ r"''/^ Bahcall-Wolf form; the normalization drops off as pc oc t~^ at late times. 
Figure [12] shows the results of a slightly more realistic calculation in a model designed to 
mimic the nearby dE galaxy M32. After reaching approximately the Bahcall-Wolf form, 
the density drops in amplitude with roughly fixed slope for r < r/i . This example suggests 
that the nuclei of galaxies like M32 or the Milky Way might have been ~ a few times 
denser in the past than they arc now, with correspondingly higher rates of stellar tidal 
disruption and stellar collisions. 

Expansion due to a central BH has been observed in a handful of studies based on 
fluid (Amaro-Seoane et al. 2004), Monte-Carlo (Shapiro & Marchant 1978; Marchant 
& Shapiro 1980; Freitag et al. 2006), Fokker-Planck (Murphy et al. 1991), and TV-body 
(Baumgardt et al. 2000) algorithms. All of these studies allowed stars to be lost into or 
destroyed by the BH; however most adopted parameters more suited to globular clusters 
than to nuclei, e.g. a constant-density core. Murphy et al. (1991) applied the isotropic, 
multi-mass Fokker-Planck equation to the evolution of nuclei containing SBHs, including 
an approximate loss term in the form of p.6p to model the scattering of low-angular- 
momentum stars into the SBH. Most of their models had what would now be considered 
unphysically high densities and the evolution was dominated by physical collisions be- 
tween stars. However in two models with lower densities, they reported observing signifi- 
cant expansion over 10^" yr; these models had initial central relaxation times of Tr < 10^ 
yr when scaled to real galaxies, similar to the relaxation times near the centers of M32 
and the Milky Way. The p ~ r"'^/'* form of the density profile near the SBH was ob- 
served to be approximately conserved during the expansion. Freitag et al. (2006) carried 
out Monte-Carlo evolutionary calculations of a suite of models containing a mass spec- 
trum, some of which were designed to mimic the Galactic center star cluster. After the 
stellar-mass BHs in their models had segregated to the center, they observed a roughly 
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Figure 12. Black-hole-driven expansion of a nucleus; this Fokker-Planck model was given pa- 
rameters such that the density at the end is similar to what is currently observed in the nucleus 
of M32, with a final infiuence radius « 3 pc. The left panel shows density profiles at constant 
time intervals after a Bahcall-Wolf cusp has been established; the right panel shows the evolution 
of the density at 0.1 pc as a function of Mace, the accumulated mass in tidally-disrupted stars. 
As scaled to M32, the final time is roughly 2 x 10^'^ yr. This plot suggests that the densities of 
coUisional nuclei like those of M32 and the Milky Way were once higher, by factors of a few, 
than at present. (From Merritt 2006b.) 



self-similar expansion. Baumgardt et al. (2004) followed core collapse in iV-body models 



with and without a massive central particle; "tidal destruction" was modelled by simply 
removing stars that came within a certain distance of the massive particle. When the 
"black hole" was present, the cluster expanded almost from the start and in an approxi- 
mately self-similar way. These important studies notwithstanding, there is a crucial need 
for more work on this problem in order to understand how the rates of processes like 
stellar tidal disruption vary over cosmological times (e.g. Milosavljevic et al. 2006). 

We thank H. Baumgardt, A. Graham and T. Lauer for permission to reproduce fig- 
ures from their published work, and A. Graham for comments on the manuscript. We 
acknowledge support from the National Science Foundation under grants AST-0420920 
and AST-0437519 and from the National Aeronautics and Space Administration under 
grants NNG04GJ48G and NNX07AH15G. 
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